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We address the issue of dephasing effects in flying polarization qubits propagating through 
optical fiber by using the method of dynamical decoupling. The control pulses are imple- 
mented with half waveplates suitably placed along the realistic lengths of the single mode 
optical fiber. The effects of the finite widths of the waveplates on the polarization rotation 
are modeled using tailored refractive index profiles inside the waveplates. We show that 
dynamical decoupling is effective in preserving the input qubit state with the fidelity close to 
one when the polarization qubit is subject to the random birefringent noise in the fiber, as 
well the rotational imperfections (flip- angle errors) due to the finite width of the waveplates. 
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I. INTRODUCTION 

Quantum information processing (QIP) has gained huge interest over the last few decades. 
This is because of the fact that it is is potentially able to solve many problems faster than the 
classical counterpart as well as provide secure communication channel. However, the inevitable 
interaction of the qubits with a noisy environment causes the loss of coherence leading to errors in 
the processing of the quantum information. This effect is known as decoherence and it limits the 
time scale over which quantum information can be retained and the distance over which it can be 
transmitted. 

Among the various strategies developed to combat the decoherence effects, we consider the 
method called dynamical decoupling (DD) [THE]. This is a relatively simple but effective technique, 
which uses sequences of external control pulses applied to the system qubits to reduce (or average 
out) the interaction of the system with the environment. A significant advantage of the DD 
techniques are that they do not require any ancilla qubits or encoding or measurement overheads. 

Most of the theoretical works on DD have hitherto considered only the case of ideal pulses. 
This means that the pulses were assumed to be instantaneous and infinitely strong in the sense 
of a (5-shaped pulse. In that case, we can ignore the effect of the noise-inducing environment 
during the application of the pulses. However, in any realistic physical implementation, this is no 
longer the case as the pulses generally have a finite duration (rp). As the effects of the applied 
pulses can be viewed as rotations of the Bloch vector on the unit sphere, the physical pulses 
result in rotational imperfections. Such imperfections typically have a cumulative effect, leading 
to progressive unaccounted phase errors in the echo signal. Such an effect leads to a considerable 
reduction of the performance of the DD sequences. 

Here, we extend the idea of dynamical decoupling with realistic pulses to long distance commu- 
nication using single photons. In optical quantum information processing, information is usually 
encoded in the polarization state of photons (which is called a polarization qubit) and photons are 
then typically routed through optical fibers or waveguides. Transmission of the photons through 
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such optical fibers can be helpful in using them as optical quantum memory [9l [10], distributed 
quantum computation [TT] and quantum cryptography [T2]. For the last few decades, quantum 
communication using propagation of the single photons through optical fibers has therefore emerged 
as a very active area of research and commercial development. The initial state of a polarization 
qubit can usually be very well prepared while it is much more challenging to preserve this state 
along the communication channel before it reaches the final detection stage. 

As the single photon propagates through the fiber, external effects such as temperature, stress, 
etc., within the fiber can randomly affect the polarization state of the photons |13j . This type 
of noise is referred to as birefringent noise as these uncontrollable fluctuating factors cause the 
birefringence, i.e. the refractive index difference An = |njy — nv\ (where nn and ny denote the 
refractive indices corresponding to the horizontal and vertical polarizations, respectively) along the 
fiber to change randomly. The effect results in the polarization state of the single photons to change 
very rapidly, making it impractical to correct for it by calibration. In practice, optical fibers used 
for communications with light can be several hundreds of kilometers long and the birefringence in 
such long optical fibers can totally destroy the information stored in the polarization qubits. 

Hence, it is crucial to protect the flying polarization qubits against such detrimental dephasing 
effects induced by the flber. Several recent studies have looked at suppression of these effects on 
the polarization qubits. Wu and Lidar [H] showed that dynamical decoupling could be applied 
for reducing quantum noise in optical fiber. Massar and Popescu presented a method to reduce 
polarization mode dispersion in optical fibers using controlled polarization rotations [15]. In our 
previous work |16| . we numerically simulated random birefringent noise along realistic fiber length 
and showed that application of DD could well preserve the input polarization qubit from such noise. 
We chose the CPMG sequence [T7] of DD in particular because this sequence has been shown to 
be robust against a variety of dephasing and control pulse errors [THHSOj; and we restricted our 
analysis to ideal pulses implemented with zero-width half-waveplates along the fiber. 

In this article, we investigate in details the issue of polarization dephasing by using finite width 
waveplates which are likely to cause some additional errors, apart from the random birefringent 
dephasing noise. Finite widths of the birefringent waveplates directly affect the phase of the photon 
transmitted through the waveplates leading to further loss of coherence. Introducing tailored re- 
fractive index profiles within the waveplates, we show that it is possible to address such detrimental 
effects with DD techniques when implemented with waveplates at the prescribed locations along 
the fiber. Estimates of the required inter-waveplate distance are provided for each of the refractive 
index profiles, and this information will be enough for an experimenter to know to successfully 
implement the faithful long-distance quantum communication channel in practice. 

This paper is organized as follows. In Sec II, we briefly introduce the basic ideas of DD (with 
ideal and realistic pulses). We discuss the nature of birefringent noise in optical fiber and how the 
DD techniques can be applied to table such noise with finite-width waveplates in Sec. III. Section 
IV contains our numerical results and comparative analysis for various tailored refractive index 
profiles. Finally, we conclude with a brief summary of the results with accompanying discussions. 

II. DYNAMICAL DECOUPLING 
A. General Idea 

Dynamical decoupling is an effective method to time-reverse the system-bath interaction by 
applying sequences of sufficiently fast and strong pulses. As a result of application of the sequence 
of pulses, the interaction of a qubit system with the environment is reduced, thus retaining the 
quantum information for longer time (or distance). The most general Hamiltonian describing the 
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FIG. 1: (Color Online) (a) Schematic representation of a cycle of general dynamical sequence with pulses of 
duration Tp, (b) CPMG sequence with two pulses in a cycle, and (c) XY-4 sequence with alternated X and 
Y waveplates. 



evolution of a system coupled to a bath can be written as 



(2.1) 



where Hs and Hb are the system and bath Hamiltonians, respectively. Since it is difficult to 
control the states of the environment, the control pulses need to act on the system, and the effect 
of these pulses are described as a refocusing of the system-environment interaction by a control 
Hamiltonian Hc{t). When the system is a qubit undergoing dephasing, we can write the interaction 
Hamiltonian Hi as 



Hi = a,0Bz, 



(2.2) 



where is the Pauli Z spin operator and Bz is a bath operator which couples to the photonic 
qubit, causing dephasing. 

Let us now consider a single cycle of a DD sequence having a period T. The evolution operator 
describing the evolution of the total system from to T (Fig. 1), in the rotating frame, can be 
written as [21] 



N 



U{T) = Uf{TN+l)llUh{Tp)Uf{n 



i=l 



where the free evolution operator is (Eq. (2.1) given by 

C/f(T) = exp[-iiJtot(T)] 
and evolution operator acting during the application of pulses is 



UhiTp) = Texp 



-t I dt {Htot + H'c{t 

'0 



(2.3) 



(2.4) 
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Here T is the standard time-ordering operator and Hamiltonian of the control pulse can be written 
as 

Hc{t) = a.fit) (2.6) 

where f{t) = {fx{t), fy{t), fz{t)) is the vector defining the shape of the pulse and a is the vector 
of the Pauli matrices. The axis of rotation due to the applied pulse at the instant t is given by 

the unit vector which implies that rotation due to the pulse can be adjusted by tuning the 

pulse shapes. The underlying principle of dynamical decoupling is to select a "pulse sequence" f{t) 
which causes the integrated time evolution of the interaction Hamiltonian to coherently average to 
zero. 



Since the evolution described in Eq. (2.3) is necessarily an unitary one, it can be written as the 
exponential of a hermitian operator H^s 

U{T) = exp[-iH,sT], (2.7) 

where i/gff can be written as a series expansion using the average Hamiltonian theory ^22j 



n=0 

An ideal DD sequence, i.e. a DD sequence with instantaneous pulses, makes H^^^ = in 
general and better performance of a DD sequence usually corresponds to progressively eliminating 
the higher order terms in such an expansion [21 [23l [21] (e.g. Magnus expansion [25])- Prominent 
examples of DD schemes are the periodic DD (PDD) |3], Carr-Purcell DD (CP) [26], Carr-Purcell- 
Meiboom-Gill (CPMG) [IT], concatenated DD (CDD) [2] and Uhrig DD (UDD) [3] 

B. Ideal and Real Pulses 

Ideally, DD pulses are assumed to be strong and instantaneous pulses applied fast enough 
compared to the internal dynamics of the environment. Under the assumption of weak coupling to 



the environment, the evolution operator in Eq. 2.5 (in the rotating frame) simplifies to 



UhiTp) = exp[-iaaep/2]. (2.9) 

Here a = x,y, z and 9p = ujpTp {ojp being the frequency of the pulses) is the rotation around the a 
axis. For ideal instantaneous pulses which implement vr rotations, the angle 9p will be vr. 

However, imperfect pulses can result in errors in the rotation axis as well as angle of rotation. 
We can write the resulting propagator as the product of the ideal pulse propagator and a rotational 
error exp[— i(Tej^ei/2] [21], due to the i-th waveplate, 

Uh{rp) = exp[-i(Je,^e,/2] exp[-m„^p/2]. (2.10) 
The modified free evolution operator in presence of the pulse errors can be written as 

Uf{Ti,Tp) = ?7/(r,)expHaeJe72]- (2.11) 



The total evolution operation from Eq. (2.3) then reads 

N 



U{T) = U'f^^^ {TN+i.Tp) n Uhm'M.Tp). (2.12) 



i=l 
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A good choice of DD sequence should make U{T) ~ I (the identity) in the presence of the 
pulse errors defined above, preserving the initial qubit against decoherence along with rotational 
imperfections. Investigating into the effects of the realistic pulses, Khodjasteh and Lidar analyzed 
the cumulative effects in pulse sequences and provided an optimum pulse interval for a fixed minimal 
pulse width Tp ^m |2Z!- Uhrig and Pasini showed the optimized performance of the DD sequences 
for considering realistic control pulses of finite duration and amplitude \28\ |29] . Composite pulse 
sequences such as BBl, CORPSE, and SCORPSE, have been shown to correct systematic pulse 
errors (which might include pulse amplitude, phase and frequency errors) [30l - [33] . Pulse shaping 
is another method that is used to counteract environmental noise effects during the finite duration 
of the real pulses [331 l35l - [37] . 



In a long distance communication channel (typically of lengths 10-1000 km), often made with 
polarization-preserving birefringent fiber, polarization qubits are likely to experience random effects 
due to changes in temperature, stress, etc. during propagation. The characteristic length scales for 
such changes may be several meters, i.e. lengths smaller than the fiber beat lengths \15\ I38j. We 
approximate the communication channel as a continuously connected fiber elements which have 
sections of constant birefringence on the order of this length scale. In the following, we consider 
the evolution of pure single photon state through such channel. 

Any arbitrary single-photon polarization state can be characterized by its polarization and 
frequency spectrum as [39] 



with \H) and \V) denoting the horizontal and vertical polarization states having arbitrary 
complex amplitudes a and /3, respectively (which satisfy |ap-|-|/3p = 1) and e{uj) being the 
complex amplitude corresponding to the frequency w. If the input photon is transmitted freely for 
a length L through the birefringent fiber, then the qubit state becomes 



Here the birefringence An{x) = \nH — nv\ causes the dephasing in the single photon polarization 
state. The phase accumulated by the qubit for the length L is given by 



III. SUPPRESSING BIREFRINGENT DEPHASING WITH DD 



A. Effect of birefringent dephasing on the polarization qubits 




(3.1) 




(3.2) 




(3.3) 



The final polarization state is obtained as 
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where the decoherence function, F{L) = J duj\e{uj)\'^ exp{iujAnL/c), gives the decay of the off- 
diagonal elements as the qubit propagates through the fiber. This analysis indicates that the 
interaction of a polarization qubit with the fiber during its propagation causes the qubit to lose 
coherence. As a result, the phase of the qubit becomes randomized and the quantum information 
stored in it is eventually lost. 



B. Characterizing the dephasing noise 

If we describe the fluctuations in the fiber to be random, stochastic fluctuation of the refractive 
index difference An can be simulated as a Gaussian-distributed zero mean random process. In this 
case, the noise is completely defined by the first order correlation function at two points xi and X2 
inside the fiber given by 



{An{xi)An{x2)) = exp [— An^/2(T^^ 



(3.5) 



Here (. . .) represents the stochastic average, i.e. with respect to the realizations of the bire- 
fringent noise and \xi — X2I is considered to be less than correlation length. (Estimates of the 
correlation lengths for a typical optical fiber are given in Ref. |3 8|.) 



The Fourier transform S{k) of the correlation function in Eq. (3.5 ) is S{k) 

is related to the decoherence function (defined as the ratio of the : 
of the density matrix) 



exp 



which 



final and initial off-diagonal terms 
In our previous work, we showed that by using the decoherence 
function formalism that it is possible to reduce the dephasing rate with the coherence preserved 
for a longer length along the fiber with a suitable DD sequence. 



C. Choice of DD to preserve the polarization qubit 

The facts that the DD pulses are implemented with waveplates in our scheme, and the axis 
of rotation is fixed by the orientation of the optical axis of the birefringent waveplates, highly 
restrict our choice of DD methods because it is technically formidable task to have precise control 
of the varied orientations of the optic axes as required in each sequence of the composite pulse 
sequences such as BBl, CORPSE, SCORPSE, or KDD (although all of them generally provide 
robust performance against pulse errors). Moreover, these sequences typically require large number 
of pulses in each cycle, which in our case of long-distance fibers affects scalability. 

In our case, we have to choose DD sequences which provide robust performance in presence 
of our dephasing model, which has randomly distributed noise with Gaussian spectral density as 



specified in Eq.( 3.5). Optimized sequence such as UDD is also not a good choice as UDD works 
best when the noise has a sharp high-frequency cut-off |20 | l40fH2] ). In the case of UDD, with single 
or multi-axis control, pulse errors generally accumulate with higher orders. A few recent studies 
also indicate that high-order UDD or concatenated DD sequences in general lose their advantage 
when the pulse intervals are strongly constrained \28\ H3] . 

CPMG and XY-4 Sequences: In order to preserve polarization states in a fiber, the CPMG 
sequence has been shown to work best in such Gaussian-distributed random birefringent noise 
\20\ m] . Another motivation for using CPMG is that this sequence is extremely robust against all 
pulse imperfections, when used with the longitudinal states, while giving marginally better results 
to preserve the transverse components of polarization [T6 t [T8H2H [3 ^ H5 | H6j . It requires vr rotations 
around a fixed axis which can be easily set by orienting the optical axis within half- waveplates. 
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For very similar reasons, it is natural to consider the XY-4 sequence (that requires alternating tt 
rotations around x and y axes), which is known to provide excellent performance in presence of pulse 
errors by preserving both the longitudinal and transverse components of polarization |2 H 1341 H7 1 148]. 
These sequences act as high-pass filters that effectively filter out the components of the Hi which 
vary slowly compared to r. In both sequences, the total evolution operator after one cycle, defined 
in Eq. (2.12), is Ut ~ I + ©(w^r^) (tc = 1/wc is the correlation time of the environment). Hence, 



the errors (due to both dephasing and pulse imperfections) resulting in the randomization in the 
phase of the polarization state coming out of the fiber at the end, can be reduced with CPMG and 
XY-4 up to the first order in ojcT for each cycle. 
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FIG. 2: (Color Online) Refractive index profiles generating the phase error: (a) Gaussian, (b) Rectangular, 
and (c) Hyperbolic tangent (as defined in the text). The distance x represents the distance within the 
waveplate in cm. 



IV. NUMERICAL SIMULATIONS AND RESULTS 

We focus only on the dephasing noise (i.e. the noise represented by the Pauli Z operator), 
and neglect the errors due to energy dissipation inside the fiber. This assumption is justified since 
we restrict our calculations for the wavelengths in the telecommunication band (around 1550 nm) 
where the optical losses are minimal. In fact, any state prepared along a direction that is different 
from the z-axis will undergo dephasing, i.e. a rotation around the z-axis, due to the birefringent 
noise along the fiber. 

Simulating birefringent noise: We model birefringent dephasing by continuously concatenating 
pieces of fiber with randomly generated lengths AL. The total propagation length thus can be 
split into segments of length AL with constant An{x). The phase difference for the i-th segment 
is equal to the sum of (27r/A)ALjAnj. These segments constitute a single phase profile associated 
with a particular realization of birefringent noise and corresponding changes in the refractive index 
difference An. Ensemble averaging over profiles gives density matrix for the output state depicting 
the random dephasing in the fiber. 

Implementing DD sequences: We investigate two DD sequences :-CPMG and XY-4 applied 
to the fiying polarization qubits propagating through optical fiber. The basic cycles of these two 
sequences are frX f2TX fr and frX Jt-Y frX frY , respectively, both with the cycle period of 4r. The 
free evolution periods correspond to the phase error accumulated during free propagation along 
the fiber for a length fr, and the X and Y rotations correspond to the vr rotations implemented 
with half waveplates, as shown in Fig. 1. We make a comparative analysis of these sequences for 
the following refractive index profiles that are used to generate the errors due to the finite widths 
of the waveplates. 
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Generating effects of the finite width: We intend to simulate the effect of finite widths of the 
waveplates when the polarization qubits are fed into the fiber. In presence of errors due to such 



waveplates, from Eq. (2.10) the total propagation operator after one cycle can be written as 



UciTp) = exp[-iae^9ej2]exp[-iaa0p/2] 



exp[ 



(4.1) 



In our model with finite-width waveplates, the error term A9p is due to the refractive index 
profile AN{x) within the waveplates (which have width Deltal) and can be written as 



AOp = (27r/A) / AN{x)dx 
Jo 



(4.2) 



We consider the following refractive index profiles for simulating realistic pulse effects: 



AN{x) 



exp 



2a^ 



■,0 < xq < Al 
1;0 < X < Al and elsewhere 
tanh[a(j: + 1) + 1] tanh[-a(x - 1) + 1]; 



(Gaussian) 
(Rectangular) 
(Hyperbolic Tangent) 



(4.3) 



For the Gaussian profile, xq and a denote the mean and the standard deviation of the refractive 
index distribution within the wave plate, and for the hyperbolic tangent profile the parameter a 
can be used to adjust the slope of the distribution. Profiles for the above refractive indices to 
simulate the finite width phase error are shown in Fig. 2. In our simulation, we considered the 
Gaussian refractive index profiles with mean xq = 2 and standard deviations up to 1.8, simulating 
thereby reasonably large flip-angle errors. For the hyperbolic tangent proflle, the parameter a, 
that determines the slope of the profile inside the waveplates, was optimized (for fixed number of 
waveplates) and the optimized value of a = 8 is used in the calculation. 




FIG. 3: (Color Online) Fidelity obtained with XY-4 waveplates in an optical fiber of length 10 km is shown 
with variation of the number of waveplates for different refractive index profiles. 



To characterize the effectiveness of our scheme, we use the fidelity J- between the input state 
iV'in) and Pout as 

^ = (V'inl/OoutlV'in), 
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FIG. 4: (Color Online) Fidelity obtained with CPMG waveplates in an optical fiber of length 10km is shown 
with variation of the number of waveplates for different refractive index profiles. 



n 



where pout = ^ Here n is the total number of randomly generated phase profiles, 

i=l 

corresponding to the propagation operator -Uj so that iV'j) = ttj|V'in) represents the simulated 
birefringent noise. Therefore, the fidelity being close to one implies that the input state is well- 
preserved against the dephasing. 

The TT rotations required for CPMG and XY-4 sequences are implemented with suitably oriented 
half waveplates, and the effects of their finite widths on the relative random phase are generated 
from the refractive index profiles of Eq. (4.3). In Fig. 4, we show how the fidelity varies and 
improves with the XY-4 sequence being applied for a realistic fiber length of 10 km, even for 
a large variation of the parameters of the random dephasing Ac/) for the chosen refractive index 
profiles. Fig. 4 illustrates the results for the CPMG sequence for the same length of the fiber. 

From these figures, we find that the while both the sequences work reasonably well to preserve 
the input polarization states for both the gaussian and hyperbolic tangent refractive index profiles, 
the rectangular refractive index profile gives the worst fidelity in both the cases. In fact, this profile 
never gives perfect fidelity with the CPMG sequence. In general, the fidelity is preserved better in 
case of XY-4 sequence (Fig. 3) than the CPMG (Fig. 4). This is due to the fact that the phase 
errors due to the finite width of the waveplates get partially cancelled due to the alternating vr 
rotations around two orthogonal optic axes (x and y) in a XY-4 sequence. It is also interesting 
to note that the fidelity in general improves with the increasing number of pulses (waveplates) in 
both cases showing the robustness of these schemes in the sense that the pulse errors tend to cancel 
each other instead of getting added up. 

The required number of waveplates to achieve a given fidelity can also be easily estimated from 
the above figures. For instance, for the hyperbolic tangent refractive index profiles, the required 
number of waveplates to achieve a 99.9% fidelity, are 840 and 860 for the CPMG and XY-4 
sequences, respectively. We also find that fewer number of waveplates are required for hyperbolic 
tangent refractive index profiles (for both the sequences) to achieve the same high fidelity, implying 
that our DD-sequences perform best when the finite width effects are simulated with this particular 
profile. 

Due to the finite widths of the waveplates, the actual angle of rotation deviates from vr and 
this constitutes the flip angle error in the polarization state of the photon. In Fig. 5, we plot the 
variations of fidelity with respect to the standard deviation of the birefringent dephasing and 
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FIG. 5: (Color Online) Contour plots of the fidelity with the variations of the standard deviations of 
the random birefringent dephasing A0 and the flip angle error for CPMG (Left) and XY-4 (Right). The 
simulations are done with fixed number of waveplates (1000) and total length of the fiber L = 10 km. 

flip angle errors for both the sequences. This shows that the input state is preserved up to fidelity 
close to one for a wide variation of the dephasing angle as well as the flip angle errors. 

V. DISCUSSION AND CONCLUSION 

For polarization qubit propagating through optical fibers, we demonstrate that the dephasing 
errors, contributed by both the fiber birefringence and the finite widths of the waveplates, could 
be suppressed by suitable dynamical decoupling methods. Regardless of the amplitude of the 
rotational error and random birefringent dephasing, our scheme provides a practical way to tackle 
them as long as the appropriate waveplate separations are maintained. As we have dealt with noises 
due to random fluctuations caused by any possible source such as temperature, stress, etc., the 
prescribed DD methods can be applied without an experimentalist having a detailed, quantitative 
knowledge of the decohering environment. To experimentally implement our proposed method 
to preserve the polarization qubits, several familiar techniques could be suitable depending on the 
range of fiber lengths one wishes to use. The waveplates may be directly incorporated into the fiber 
during the manufacturing process. Other methods include writing a Bragg transmission grating 
periodically into the fiber [50], or twisting the fiber in controlled ways causing suitable mechanical 
stress \13\- 

We have successfully showed that it is possible to combat the random birefringent noise in 
an optical fiber with DD waveplates which have finite widths. Our approach provides a practical 
approach to minimize errors due to both random dephasing noise in the fiber and the rotational 
errors due to implementation of real finite-width waveplates. This will be helpful to improve the 
range of quantum communication channels without requiring expensive resources such as ancilla 
or measurements and hence it leads to immediate commercial applications for quantum telecom- 
munication with light. The control overhead in the proposed application of the DD sequences 
being reasonably small, our scheme will reduce the dephasing error while implementing a scalable 
quantum computing scheme with photonic qubits. 
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